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ABSTRACT 

Thanks to the recent discovery by Fermi of about fifty new gamma-ray pulsars, it 
becomes possible to look for statistical properties of their pulsed high-energy emis- 
sion, especially their light-curves and phase-resolved spectra. These pulsars emit by 
definition mostly gamma-ray photons but some of them are also detected in the radio 
band. For those seen in these two extreme energies, the relation between time lag of 
radio/gamma-ray pulses and gamma-ray peak separation, in case both high-energy 
pulses are seen, helps to put some constrain on the magnetospheric emission mecha- 
nisms and location. This idea is analyzed in detail in this paper, assuming a polar cap 
model for the radio pulses and the striped wind geometry for the pulsed high-energy 
counterpart. 

Combining the time-dependent emissivity in the wind, supposed to be inverse 
Compton radiation, with a simple polar cap emission model along and around the 
magnetic axis, we compute the radio and gamma-ray light-curves, summarizing the 
results in several phase plots. The phase lag as well as the gamma-ray peak separation 
dependence on the pulsar inclination angle and on the viewing angle are studied. Using 
the gamma-ray pulsar catalog compiled from the Fermi data, we are able to predict the 
radio lag/peak separation relation and compare it with available observations taken 
from this catalog. 

This simple geometric model, combining polar cap and striped wind radiation 
performs satisfactorily at explaining the observed radio/gamma-ray correlation. This 
supports the idea of distinct emission locations for the radio and gamma-ray radia- 
tion. Nevertheless, time retardation effects like curved-space time and magnetic field 
lines winding up close to the neutron star can lead to discrepancy between our pre- 
dicted time lag and a more realistic relation as deduced from the gamma-ray catalog. 
Moreover, no accurate polar cap description being at hand so far, large uncertainties 
remains on the altitude and geometry of the radio emission. 

Key words: Pulsars: general - Radiation mechanisms: non-thermal - Gamma rays: 
observations - Gamma rays: theory - Stars: winds, outflows 



1 INTRODUCTION 

The new catalog on ga mma-ray pulsars o btained by the 
Fermi-LAT instrument (|Abdo et alj [2010a) increased the 
number of gamma-ray pulsars from seven to about fifty. 
Since then, new pulsars are discovered regularly. This al- 
lows for the first time a reasonable statistical analysis on 
the high-energy emission properties of these objects like 
spectral shape, cut-off energy, and comparison between ra- 
dio and gamma-ray radiation if available. For some high 
luminosity pulsars, this analysis has b een completed by a 
phase-resolved study like for the Crab (|Abdo et alj|2010bl ). 
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Vela (|Abdo et alj|2009h and Geminga (|Abdo et al.ll2010d ). 
Radio pulses and gamma-ray photons are expected to be 
produced in different emission sites, probably close to the 
neutron st ar surface for the former, des cribed by a polar 
cap model ( Radhakrishn an" fc Cookell 19691 ). and in the vicin- 
ity o f the light-cylinder for the latter, explained by outer 
gaps IrCheng et aLll 19861 ) or alternatively by the striped wind 
IjPetril 20091 ). Polarization properties at multi- wavelength 
would cer tainly help to con s train the geometry of the sites o f 
emission (|Dvks et alj|2004 iPetri fc Kirkll2005l : |Petrill2009l ). 
Gamma-ray light-curves alone can already give go od insight 
into the magnetosphere (|Romani fc WattersH2010l ). 

The new sample of Fermi-LAT gamma-ray pulsars in- 
creased interest into modeling of gamma-ray emission. De- 
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tection of several millisecond gamma-ray pul sars was not 
expec ted and came as a real surprise. Thus, IVenter et al.l 
(2009) focused special attention to this class of millisecond 
pulsars to probe the geometry of the emission regions, taking 
into account relativistic effects. 

Very recently, gamma-ray light-curves have 
been computed for the simple vacuum dipole model 
(|Bai fc Spitkovskvll2010bh or better for a realistic magneto- 
spheric model based on 3D MHD simulations of the near 
pulsar magnetosphere. This requires some post-processing 
prescription about the e mission location and me chanism 
within the magnetosphere l|Bai fc Spitkovskvl2010al ). In this 
model, gamma rays are expected close to the light-cylinder. 

An alternative site for the production of p ulsed radi- 
ation has been investigated a few years ago by iKirk et ail 
(2002). This model is base d on the s tripe d pulsa r wind , 
originally introduced bv lCoronitil (|l990l ) and lMichell |l994l ). 
Emission from the striped wind originates outside the light 
cylinder and relativistic beaming effects are responsible for 
the phase coherence of this radiation. It has already been 
shown that this model can satisfacto rily fit the optical po- 
larization data from the Crab pulsar (jPetri fc Kirkll2005l ) as 
well as th e phase-reso lved high-energy spectral variability of 
Geminga (|Petrill2009h . 

The aim of this work is to compl ete the atlas o f light - 
curves and phase plots performed by IWatters et al l (|2009r i 
on the ground of the polar cap, the slot gap and the outer 
gap models. It furnishes a extended benchmark to test dif- 
ferent scenarios for the physical processes occurring in the 
pulsar magnetosphere. We use an explicit asymptotic so- 
lution for the large-scale m agnetic field stru cture related to 
the oblique split monopole (jBogova lov 1999) and responsible 
for the gamma-ray light-curves combined with a simple po- 
lar cap geometry for the radio counterpart. Details are given 
in Sec. [2] We then compute the properties of the radio and 
high-energy light curves for the pulsed emission and com- 
pare our results with several gamma-ray pulsars extracted 
from the Fermi catalog. Relevant results are discussed in 
Sec. [3] before concluding this study. 



2 THE MODEL 

The model employed to compute the gamma-ray pulse shape 
emanating from the striped wind is briefly reexamined in 
this section. The geometrical configuration and emitting 
particle dist ribution func tions follows the same lines as those 
described in iPetril (|2009l ). The magnetized neutron star ro- 
tates at an angular speed of fi* (period P = 2-k/Q,*) di- 
rected along the (Oz)-axis, i.e. the rotation axis is given 
by f2* = f2, e z . We use a Cartesian coordinate system with 
coordinates (x,y,z) and orthonormal basis (e x , e y , e z ). The 
stellar magnetic moment is denoted by it is assumed to 
be dipolar and makes an angle x with respect to the rotation 
axis such that 

/!» = fi t [sinx (cos(f2* t) e x + sin(f2» t) e y ) + cosx e z ]. (1) 

This angle is therefore defined by cosx — ' Cz/a 1 *- The 
inclination of the line of sight with respect to the rotational 
axis, and defined by the unit vector n, is denoted by It 
lies on the (Oyz) plane, thus 



Thus we have cos£ = n ■ e z . Moreover, the wind expands 
radially outwards at a constant velocity V close to the speed 
of light denoted by c. 

Because we are not interested in the phase-resolved 
spectra, but only in the light-curves above 100 MeV, our 
model only involves geometrical properties related to the 
magnetic field structure and viewing angle. The particle dis- 
tribution function does not play any role except for its den- 
sity number. Dynamical properties (such that energy dis- 
tribution) related to the emitting particles are unimportant 
in this case. In order to compute the light curves, we use a 
simple expression for the emissivity of the wind related to 
inverse Compton scattering. This is explained in the next 
paragraphs. 



2.1 Magnetic field structure 

The geometrical structure of the win d is based on the a symp- 
totic magnetic field solution given bv lBogovalovl l|l999h . Out- 
side the light cylinder, the magnetic structure is replaced by 
two magnetic monopoles with equal and opposite intensity. 
The current sheet sustaining the magnetic polarity rever- 
sal arising in this solution, expressed in spherical coordi- 
nates (r, if) is defined by 



r s (#, <p, t) — f3 v rL ± arccos(— cot 1? cot x) H ip + 2ln 

(3) 

where /3 V = V/c, — c/S7* is the radius of the light cylin- 
der, t is the time as measured by a distant observer at rest, 
and I an integer. Because of the ideal MHD assumption, 
this surface is frozen into the plasma and therefore also 
moves radially outwards at a speed V. Strictly speaking, 
the current sheets are infinitely thin and the pulse width 
would be inversely pro portional to the wind Lorentz fac- 
tor r v = (1 - Ply 1 ' 2 , (|Kirk et al.ll2009T l. Here, as already 
done for the st udy of the synchro tron polarization of the 
pulse d emission l|Petri fc Kirkll2005l ) and for Geminga (jPetril 
2009) we release this restrictive and unphysical prescription. 
Indeed, the current sheets are assumed to have a given thick- 
ness, parameterized by the quantity A v (see Eq. ([6]) below 
for an explicit expression of the thickness). Moreover, inside 
the sheets, the particle number density is very high while the 
magnetic field is weak. In whole space, the magnetic field is 
purely toroidal and given by 



„ Rl 

r 



(4) 



The strength of the magnetic field at the light-cylind er is 
denoted by Bl • In the original work of lBogovalovl (|l999T ) , the 
function r\ v is related to the Heaviside unit step function and 
can only have two values ±1, leading to the discontinuity in 
magnetic field. In order to make the transition more smooth, 
we redefine the function r\ v by 



tanh(A v ip) 



(5) 



tp = cos i? cos x + sm $ sm X cos \<P ~ ^* ~ y 7 ) J 

With these formulae, the physical length of the transition 
layer has a thickness of the order of 



sin £ e y + cos £ e z 



(2) 



2 7r/3 v r L /A ¥ 



(6) 
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2.2 Particle density number 

The aim of this paper is to show the behavior of the pulsed 
high-energy light-curves emanating from the striped wind 
flow for different magnetic plus emitting particle configura- 
tions. We will not perform a detailed study of the phase- 
resolved spectral variability. Thus for this purpose, it is suf- 
ficient to fix the particle density number without specifying 
the distribution in momentum space. We will assume a e 
pair mixture flowing dominantly within the current sheets 
defined by Eq. |(3]). We thus adopt the following expression 
for the total particle density number 



K e (r,t) = 



(N-N )sech 2 (A v iP) + N 



(7) 



No sets the minimum density in the stripes, between the 
current sheets, whereas N de fi nes th e highest density inside 
the sheets. We refer to IPetril (|2009l ) for more details about 
the justification of this choice. 



2.3 Inverse Compton emissivity in the wind 

These particles will be visible through their inverse Comp- 
ton radiation. The corresponding total emissivity is denoted 
by j°c B - Strictly speaking, this radiated power should de- 
pend on space location f, time of observation t as well as on 
the energy of the detected photons. However, recall that a 
spectral study is out of the scope of this work, so the energy 
dependence can be dropped. We only retain variation with 

f and t. 

As done in IPetril (|2009l l. the inverse Compton light 
curves are obtained by integrating this emissivity over the 
whole wind region. This wind is assumed to extend from a 
radius ro to an outer radius r s which can be interpreted as 
the location of the termination shock. Therefore, the inverse 
Compton radiation at a fixed observer time t is given by 



7° bs (i) 



ii° bs (r, tret) r 2 sin 'ddrd'&dif 



(8) 



r J J 



The retarded time is expressed as t re t = t — | \Ro — r\ \/c w t — 
Ro/c+n-r/c. Eq.([8| is integrated numerically. We compute 
the inverse Compton intensity for several geometries. We are 
therefore able to predict the phase resolved pulse shape, for 
any inclination of the line of sight and any obliquity of the 
pulsar. Results and applications to some 7-ray pulsars are 
discussed in the next section. 



2.4 Polar cap emissivity 

The polar caps or their close vicinity to an altitude up to a 
few stellar radii from the stellar crust are the favorite sites 
to explain the coherent pulsed radio emission. For our ge- 
ometric model, both the polar caps radiate most efficiently 
at their center which means along the magnetic axis. When 
moving away from the poles staying on the stellar surface, 
the radio intensity should decrease. We therefore choose a 
smooth gaussian decreasing with distance from the magnetic 
north or south pole, almost complete extinction occurring 
outside the polar caps. Their size is determined by the loca- 
tion of the foot of the last open magnetic field lines, on the 
stellar surface. For simplicity, neither bending of magnetic 
field lines nor (general-)relativistic effects nor plasma effects 



are included. The distortions implied by these effects are not 
included in our description. Moreover, using the expression 
for the aligned rotator, the size of one polar cap, assumed to 
have a circular shape, is obtained from its angular opening 
angle # po leading to a radius R pc such that 




(9) 
(10) 



The maximal intensity arises when line of sight and mag- 
netic moment are aligned or counter-aligned. For a sharp 
transition between on and off states, emission is expected 
only when the angle formed by /?* and ft is less than the 
angular opening angle, thus (/2»,n) ^ # pc . In other words, 
a polar cap (either north or south) is only seen if 



cos (/2* , ft) = 



fj,* 



ft, 



^ cos tf pc = J 1 



(11) 



This is translated mathematically for our gaussian by intro- 
ducing the intermediate function & pc defined as 



*pc(*) = 



= cos £ cos x + sin C sin X sin(f2„ t) (12) 



Therefore one polar cap is visible if 

KM > 1 - ^ (13) 

Expressed as a polar cap emissivity e pc with gaussian similar 
shape, we write it 



<(*o 



(14) 



with a a positive constant controlling the extension of sig- 
nificant polar cap emission. e pc is maximum for /2„ and ft 
colinear with value close to unity and tends to zero for large 
distance from the magnetic poles. We stress that this de- 
scription is purely phenomenological, to be included in the 
phase plot for gamma-ray light-curves. 



3 RESULTS 

Our simple geometrical model allows to perform some ana- 
lytical calculations to deduce the relative lag in arrival time 
between radio and gamma-ray photons and the gamma-ray 
peak separation, if both, radio emission and double peak 
structure are detectable. In the following paragraphs, we 
perform a detailed analysis of the light-curve properties and 
correlation with radio wave-bands depending on \ an d C- We 
then close the discussion with application to some gamma- 
ray pulsars from the Fermi first year catalog. 



3.1 The high-energy peak separation 

In the striped wind model, the light-curve for the pulsed 
emission shows usually (in a favorable configuration) a dou- 
ble peak structure. In this case, there exist a simple relation 
between this double gamma-ray pulse separation, denoted 
by A, and the geometry of the system, assuming a rotating 
dipole with magnetic obliquity \ an d une °f sight inclina- 
tion (. Indeed, to understand the physical mechanism, let 
us assume that the wind radiates mostly when crossing a 
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spherical shell located at a radius _Ro- Due to the strong rel- 
ativistic beaming effect for ultra relativistic expansion, the 
observer only sees light travelling towards him, i.e. in the di- 
rection ($ = (p = 7r/2). For particles located in the current 
sheets, this corresponds to an observer time t a satisfying 



Ro = r B (d,Lp,t< 



(15) 



± arccosf — cot C cot y) + — — + 2 1 n 
r L 2 



where I £ Z is a natural integer. Solving for this time t a at 
which the photon leaves the current sheet, we find 



Ro arccos(— cot£ cot x) 



7-" 



(16) 



c/3 v ft* 

The ± sign labels the arm of the double spiral structure (seen 
in the equatorial plane) responsible for the radiation. Recall 
that the double peak light-curves follow from this double 
spiral shape. The time separation between two consecutive 
pulses becomes therefore 

arccos(— cot £ cot x) 



At — r 



t: = 2 ■ 



(17) 



Note that this is not necessarily the high-energy peak sep- 
aration because it can in principle vary between zero and a 
full period P. To be more precise, we have to compare one 
specified pulse with its two neighbors in time, the one com- 
ing earlier to the observer and the other later. Doing so we 
get 



At 



peak 



minflC-*. |,|C - 
min(|At cons |, |At cons 



(t7+P)\) 
~P\) 



(18) 
(19) 



Normalizing to the period of the pulsar P, the phase sepa- 
ration reads 



A = 



At 



,„,,]. . /arccos(-cotC cotx) 
= mm 2 — . 



P v 
arccos(— cot £ cot x) 



From the definition of the arccos function, we get 

s(— cot C cot v) ■£ , > , ^ n 

— ^ if cot Q cot x ^ 

if cot £ cot x ^ 



arccosf — cot £ cot x) 



(20) 
(21) 

(22) 



Finally, the relation between obliquity x, line of sight in- 
clination angle £ and peak separation A can be rearranged 
into 



cos(-7r A) = 



- cot £ cot x if cot £ cot X ^ 



(23) 



cot £ cot x if cot £ cot x ^ 

To summarize we find 

cos(7r A) = | cot C cot x! (24) 

This relation is shown in Fig. [1] where cos ( is plotted versus 
cos x for constant value of the separation A starting from 
zero to 0.5 with a step 0.02. 

A few special cases are worth to mention. First, zero 
separation or overlapping of both pulses, corresponding 
to A = 0, implies £ = n/2 — \- For the chosen variables 
in Fig. [T] it corresponds to the quarter circle of radius unity 
centered at the origin. Second, a separation of half a period 
A = 0.5, implies a line of sight contained in the equato- 
rial plane of the pulsar, £ = 7r/2. It is independent of the 
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Figure 1. Peak separation A between both gamma-ray pulses 
purely determined by the obliquity of the pulsar \ and the incli- 
nation angle of the observer line of sight £. Each curve represents 
a constant value of A ranging from to 0.5 by 0.02 steps. No sep- 
aration occurs along the circle of radius unity (pulse overlapping) 
whereas maximum separation of 0.5 is reached for cosf = 0. The 
diagonal f = \ describes the visibility of one magnetic pole. 



obliquity. This corresponds to the straight horizontal line 
along the x-axis in Fig . [1] The relation Eq. (|24|l was already 
noticed by Kirk (2005), see figure therein. For symmetry rea- 
sons, we will assume that \ e [0, 7r/2], therefore cotx ^ 0. 



3.2 The radio lag 5-peak separation A relation 

Next, we show that the combination polar cap/striped wind 
allows to derive a simple analytical relation between the lag 
of radio vs gamma-ray arrival time <5 and the high-energy 
peak separation A. Because the two sites of emission are 
very distinct, inside and outside the light cylinder, the lag 
is interpreted as a time of flight delay between the region of 
radio radiation and the one for gamma-ray photon produc- 
tion. 

The choice of origin of time and line of sight located 
in the plane {Oyz) (therefore ip = ^) implies that the radio 
photon emitted towards the observer will happen at periodic 
dates with period P, for one pole, given by 



t 



TV , 2 7T 

— — + k — 

2f2* ft* 



+ kP 



(25) 



where k £ Z is a natural integer. The superscript n denotes 
bundles of photons emanating from the same polar cap, let 
us say the north magnetic pole. The same applies to the 
south cap but due to the symmetry, the phase difference 
corresponds to half a period, P/2, thus 



P 3P , n 



(26) 



The observer remains at a distance D from the center of 
the neutron star. The polar cap photon needs therefore a 
time At = (D — R*)/c to reach him after emission from the 
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pulsar. Thus the detection time is the sum 

D — R, 
ti — t„ c . H 



Radio time lag 



(27) 



The gamma-ray photon from the current sheets arrives ear- 
lier because closer to the observer, he needs a shorter flying 
time given by At = (D — Rq)/c leading to a time of detection 
given by 



± D 
t 2 =tt + — 



Ro 



(28) 



As a consequence, the difference in arrival time between both 
kind of photons, focusing on the radio emission coming out 
of the north pole (visible for < £ < 7r/2), is 



fa 



t2 — tl = t s 



R» - Ro 



1 — /3 V Ro _j_ arccos(— cot £ cot x) , R* 
/? v c Q* c 



(29) 
(k + l)P 



By normalizing this time to the period of the pulsar P and 
using Eq. (|22|) with cot £ cot x ~& 0, we arrive at 



ta 

P 



Ro 



/3y 27rr L 



± ■ 



1 - A 



2 7rr L 



(fc + Z) 



(30) 
(31) 



The same procedure applied to the south pole with n/2 < 
(" < 7r implies cot £ cot x ^ 0, leading to the same conclusion 
because of symmetry. 

Recall that in order to observe pulsed emission from the 
radially expanding wind, the flow has to be relativistic with 
T v 3> 1 and Ro must lie close enough to the light-cylinder 
such that the underlying condition is fulfilled 



Ro <2TYT 2 v r L 
Thus the upper limit to the first term in Eq. (I31|l is 



I3 V Ro 



1 -Ro 
/3 V 27rr L " 2T? 27rr L 



(32) 



(33) 



Moreover, for all the pulsars in the catalog, even for the 
millisecond pulsars, the light-cylinder radius is much larger 
than the stellar radius, a > i?,. To a good approxima- 
tion, we can neglect the term 2 ^r L in Eq. (|3ip . the neutron 
star is located well within its light-cylinder. The remaining 
expression for a positive lag is 



1 - A 



(34) 



This has to be comp ared with the 8~ A diagram published in 
IXbd o et al ] (|2010ah . The data points as well as the relation 
Eq. (|34[) are shown in Fig. [2] All the points but one lie below 
the line given by Eq. (|34|l . The spread in radio lag along 
this line can be understood as a fluctuation in the precise 
location of the most significant emission in the wind. 

Finally, we give an expression for the linear fit to the 
sample of pulsars showing double peaked gamma-ray light- 
curves simultaneously with radio emission by 



Am ~ 0.59 - 1.06(5 
shown in blue line in Fig. [2] 



(35) 
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Figure 2. Gamma-ray peak separation A vs radio time lag S as 
observed by Fermi-LAT. The blue curve is a linear fit and the red 
straight line the upper limit as predicted by Eq. 1341 1 . 



3.3 Other retardation effects 

Including general-relativistic effects do not change much the 
time lag between radio and gamma-ray photons. Indeed us- 
ing the Schwarzschild geometry, for a photon with radial 
motion, the time required to reach a distant observer lo- 
cated at a distance D, starting on the stellar surface is 



cAtp C = D - R, + R s In 



D-Rs 
R» — R s 



(36) 



R s = 2 G M/ c 2 being the Schwarzschild radius. A similar 
expression applies for gamma-rays 



c AU = D - Ro + R s In 



D-Rs 
Ro — Rs 



(37) 



Therefore the time of flight difference compared to Newto- 
nian gravity is 



a j A , A , R* — Ro , Rs , ( R* — Rs 

Ai G R = At 7 - At pc = 1 In — — 

c c V Ro - Rs 



(38) 



Normalized to the pulsar period we get an increase by an 
amount of 



c Rs , I R* — R s 
ogr = hi 



2 7rr L 



Ro — Rs 



(39) 



For neutron star parameters, this remains negligible com- 
pared to unity because R a <C n,. Typical values are <5gr ~ 
0.01. 



3.4 Phase-inclination diagram 

Next, we show the phase-inclination diagram for different in- 
clinations of the line of sight £ as well as different magnetic 
obliquities x °f the pulsar. This is done for the high-energy 
photons as well as for the radio emission coming from the 
polar cap, see Fig. [3] In order to point out the relative phase 
shift between both emission processes, we normalized inde- 
pendently their respective maximum intensity to unity. The 
continuous rotated S-shaped structure is a characteristic of 
the striped wind whereas the two spots correspond to a map- 
ping of the two polar caps. Both structures are symmetric 
with respect to the equatorial plane, i.e. compared to an 
inclination £ = 90°, as expected from the geometry. 
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Figure 3. Phase plot of the pulsed gamma-ray and radio emission components for a full period of the pulsar (phase G [0,1]) for an 
inclination of the line of sight f between and 180 degrees and magnetic obliquities \ = 18, 36, 54 and 72 degrees, from top left to 
bottom right. Note the colour-coded range shown in the legends of each phase plot: white colour corresponds to the faintest phase (off 
pulse) with intensity around 0.11 whereas red colour corresponds to the brightest luminosity i.e. 1 (in normalized unities). 



The observer will see different kind of light-curves de- 
pending on his viewing angle. First, radio photons are ob- 
servable only when the line of sight intersects at least one 
polar cap, in the vicinity of the magnetic poles, that is when 
£ w X ± $ pc . Second, to observe the gamma-ray pulsation, 
this same line of sight should intersect the current sheets, 
or in other words tt/2 — \ ^ £ ^ 7r/2 + X- Conversely, no 
gamma-ray pulses are seen when ( ^ n/2~x or C ^ n /2 + X- 
Thus a necessary condition to detect simultaneously radio 
and gamma-ray photons is x ^ 7r /4 and C ~ X- F° r those 
particular pulsars, possessing a double gamma-ray peak and 
radio emission, we are able to find the geometrical param- 
eters from the peak separation, see Fig. [1] Third, a single 
pulse in the high-energy light-curve occurs when the line of 



sight passes just through the edge of the striped part of the 
wind, this can be interpreted as the special case £ ~ n/2 — x- 
In the very special case of one gamma-ray pulse and observ- 
able radio emission, this leads immediately to the geomet- 
rical parameters ( ~ x ~ t/4 = 45°. This seems to be the 
case for PSRJ0437-4715 and PSRJ2229+6114. 

The phase plots shown in Fig. [3] have been computed 
for different inclination angles £ and obliquities x but with 
constant current sheet thickness, parameterized by A v , and 
constant particle density contrast, which means constant 
ratio N/No. More precisely, we adopted A v — 10 and 
N/No = 10. This last value explain the factor around 10 in 
intensity between off-pulse and on-pulse phase. Indeed, by 
decreasing the fluctuation in particle density between the 
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hot unmagnetized and cold magnetized part of the wind, as 
done in Fig. [4] we recognize a similar trend in the intensity 
diagram. If the ratio N/No = 5, the luminosity varies be- 
tween 0.22 ~ 1/5 and 1 in normalized unities whereas for 
N/N = 2, it lies between 0.54 w 1/2 and 1. 

Moreover, the current sheet thickness directly impacts 
on the duty cycle of the light-curve, or in other words, on 
the width of the gamma-ray pulses. This has been check by 
changing the parameter A v to 5 or 2 instead of the previous 
value of 10. Results are shown in Fig. El 

To better assess the differences, we summarize all the 
light-curves for x = 72 degrees and f = 90 degrees in two 
plots as depicted in Fig|6] 

Finally, there exist two limiting cases. First the aligned 
rotator showing no pulsed emission. Second the perpendicu- 
lar rotator emits pulsed high-energy radiation over the whole 
sky and both its polar caps are visible. A direct consequence 
is a double peak structure in both radio and gamma-ray 
light-curves, with a separation of A = 0.5. PSRJ0030+0451 
is a typical example. 




20 40 60 80 



in degree 

Figure 7. Flux correction factor /n for the striped wind model vs 
obliquity \ an d inclination of Earth line of sight C,e- Two different 
regions with roughly constant factor are clearly distinguishable 
along ( E = tt/2 - X - 



3.5 Flux correction factor 

The photon flux Fobs measured by an observer located on 
Earth is biased due to anisotropic emission from the wind 
depending on the viewing angle f . This fact is clear from the 
aforementioned phase-inclination plots shown in the previ- 
ous paragraphs. The observed gamma-ray flux has to be 
corrected to obtain the true gamma-ray luminosity by in- 
troducing a correction factor fa defined by 



Lj =4 7r/ n F obs £> 2 



(40) 



Here D is the distance of the pulsar to the observer and F b s 
the obs erved flux. As in th e polar cap and outer/slot gap 
models dWatters et al.ll2009l ). the correction implied here by 
a relativistic beaming effect is given by 



fo(x,Cn 



Jo Jo *_ F i (X> f) sin j dl f 
2 Jo" F j{x,CE,<p)d<p 



(41) 



For the striped wind model, this correction factor is shown 
in fig. with the full dependence on obliquity \ and in- 
clination of Earth line of sight C,e- We can approximately 
separate the correction factor into two regions of constant 
value. In the first region, for an obliquity x < t/2 — Ce, the 
correction is close to be uniform and equal roughly to 1.90. 
This case corresponds to a line of sight not crossing the cur- 
rent sheets in the wind, there is almost no pulsed emission 
visible. In the second region, where \ > ""/2 — Cei the cor- 
rection is also almost uniform and equal roughly to 0.4. This 
case corresponds to a line of sight intersecting the current 
sheets, leading to pulsed emission. This behavior is expected 
from the definition of the correction factor Eq. (I41|l . Indeed, 
for a fixed obliquity \i t fle numerator is a constant whereas 
the denominator depends on the line of sight towards the 
Earth. On one side, in the first region \ > t/2 — £e, the 
phase-averaged emission is faint and weakly pulsed. It fol- 
lows a small denominator therefore a large correction factor. 
On the other side, in the second region, the situation is op- 
posite, the emission during the pulses is strong and so the 
denominator of Eq. (|41[) larger, finally the correction factor 
is weakest. 



Note that these bounding values depend on the sheet 
thickness parameterized by A v as well as on the particle 
density contrast parameterized by iVo and N. 



3.6 Light-curve fitting and gamma-ray luminosity 

We conclude our study by fitting light-curves of a small 
sample of pulsars. The only relevant free parameters in our 
model are the geometry of the wind, the particle density 
number and the size of the current sheets. In this last section, 
we specialize our results to some Fermi detected gamma- 
ray pulsars and show the best parameters fitting their high- 
energy light-curves above 100 MeV. Therefore, the knowl- 
edge of the viewing angle and the obliquity allows an estima- 
tion of the flux correction via the beaming factor. Eventually 
a true gamma-ray luminosity versus spin-down luminosity 
can be plotted. 

We start with an estimate of the peak separation when 
one radio pulse is detected. In that case, £ m x as explained 
above and the knowledge of A immediately implies a solu- 
tion for £. This has been done for several pulsars and listed 
in Tab. [T] In all the results shown below, for simplicity, we 
took a constant Lorentz factor of the wind equal to T v = 10. 

Let us have a deeper look on a representative sample of 
Fermi-detected gamma-ray pulsars. 



3.6.1 PSR J0030+0451 

PSR J0030+0451 is a millisecond pulsar, P = 4.87 ms, show- 
ing a double pulse structure in the radio band. This may 
suggest that both its magnetic poles are visible, or in other 
words, it is almost a perpendicular rotator with \ ~ 90°. 
Nevertheless, the maximal intensity of both radio pulses are 
sensibly different, this is interpreted as the line of sight pass- 
ing closer to one magnetic pole than to the other. An alter- 
native would be to explain it by a process occurring in the 
vicinity of the polar caps with variable efficiency. However, 
polar cap emission is not the main purpose of this work 
so we simply assume identical shapes for both radio-pulses. 
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Figure 4. Phase plot of the pulsed gamma-ray and radio emission components for a full period of the pulsar (phase G [0,1]) for an 
inclination of the line of sight £ between and 180 degrees and a magnetic obliquity \ = 72 degrees. Only the ratio N/Nq changes. 
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Figure 5. Phase plot of the pulsed gamma-ray and radio emission components for a full period of the pulsar (phase G [0,1]) for an 
inclination of the line of sight £ between and 180 degrees and a magnetic obliquity x = 72 degrees. Only the parameter Aip varies. 



Therefore, an obliquity close to 90 degrees but less matches 
the right geometry. Setting \ ~ 67° satisfactorily agrees 
with the radio light-curve, see Fig. O on the top left plot. 
To have the right gamma-ray peak separation, we have to 
adopt £ ~ 85°. Moreover, the radio-pulses are very broad, 
each of them having a duty cycle of roughly 0.2 in phase. 
Therefore, to reconcile our model with data, we have to ex- 
tend the polar cap region to a sizeable fraction of the whole 
neutron star surface, Fig. [5] top left plot. Emission starts 
right after the light- cylinder radius. There is still a excess of 
0.1 in the phase delay compared to observation. This has to 
be explained by some other retardation effects of the radio 



pulse such as the strong gravitational field regime (which we 
have shown to be negligible) or by magnetic field bending 
due to charges flowing within the magnetosphere and dis- 
turbing the closed field lines structure taking to be an exact 
dipole. 



3.6.2 PSR J0218+4232 

PSR J0218+4232 is another millisecond pulsar, P = 
2.32 ms, showing a less clear double pulse structure in 
gamma-rays. Its radio-pulse shape looks much more compli- 
cated with something like conal and core component, over 
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Figure 6. Light-curves of the pulsed gamma-ray emission component, for a full period of the pulsar (phase G [0, 1]) and for an inclination 
of the line of sight £ = 90 degrees and a magnetic obliquity x = 72 degrees. On the left picture, for different parameters N/No with 
Atp = 10 and on the right picture for different A v with N/Nq = 10. 
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Table 1. Geometry of the pulsar wind for 8 gamma-ray pulsars. The angle of inclination and obliquity are given in degrees. The Lorcntz 
factor of the wind is the same for the whole sample and taken to be T v = 10. 



a wide range of the pulsar period. The best fit is obtain for 
an obliquity \ = 57° and an inclination of C, = 75°. The 
gamma-ray light curve adjusts well to Fermi data. Moreover 
the gamma-ray pulse time lag compared to the middle of the 
radio pulse matches precisely the measurements. Note that 
the gamma-ray off pulse emission remains at an appreciable 
level over the full period, see Fig. [8] on the top right plot. 

3.6.3 PSR J 0437 -47 15 

PSR J0437-4715 is a special case of millisecond pulsar, P — 
5.76 ms, showing only one gamma-ray pulse combined with a 
sharp radio pulse. As explained in a previous discussion, the 
unique geometry allowing such a behavior needs an obliquity 
almost equal to the inclination of the line of sight, namely 
( w x « 45°. Indeed, we take \ = 45° and £ = 40° and 
arrive at the light curve presented in Fig. second row, 
left plot. For one gamma-ray pulse, our model predicts a 
delay of 0.5 in phase between radio and gamma-ray pulsars, 
in relative good agreement with Fermi measuring 0.43. The 
overestimate is about 0.07. 



3.6.4 PSR J1124-5916 

PSR J1124-5916 shows a sharp double gamma-ray peak with 
significant off-pulse emission in combination with a single 
radio peak. Adopting the parameters x — 80° and £ = 85°, 
our model is compared with observation in Fig. [8] second 



row, right plot. The expected time lag 0.255 is not very 
different from the measured 0.23. 

3.6.5 PSR J 2021 +36 51 

PSR J2021+3651 gamma-ray light-curve shows a narrow 
double pulse structure and a single radio pulse. The peak 
separation is close to 0.5. Here again, for the best fit pa- 
rameters x = 73° and £ = 65°, the predicted lag of 0.265 
overshoots the observation by roughly 0.1, Fig. [8] third row, 
left plot. 

3.6.6 PSR J2032+4127 

PSR J2032+4127 is an example of half-a-period peak separa- 
tion implying a line of sight inclination of £ < 90°. Moreover, 
because we see radio pulse it should be almost a perpendicu- 
lar rotator. However, this would permit to detect both radio 
component which is not the case, so we must conclude that 
the inclination is slightly different from an orthogonal ro- 
tator with a small polar cap reproducing a width of 0.1 in 
the radio pulse, Fig. [8] third row, right plot. We found best 
agreement for x = 89° and £ = 70°. 

3.6.7 PSR J2043+2740 

It seems that PSR J2043+2740 shows evidence for three 
gamma ray pulses, Fig. [8] bottom left plot. The two extreme 
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ones are significant whereas the middle one is less significant. 
We tried to fit the light curve according to the two well 
separated pulses. The time delay is overestimated by 0.12. 

3.6.8 PSR J2229+6114 

Finally, PSR J2229+6114 is another example of single pulse 
gamma-ray pulsar. The best fit parameters correspond to 
the same values as for PSR J0437-4715. Now, the time lag 
is even better, again the model predicts 0.5 compared to the 
observed 0.49, Fig. [SI bottom right plot. 

Details of the parameters used to fit the data are sum- 
marized in Tab. [T] 

This ends the sample of fitted gamma-ray pulsars. We 
demonstrated that our striped wind/polar cap model im- 
plies a well-defined relation between gamma-rays and radio 
pulses. Single and double pulse light-curves are expected. 
For some pulsars, the model does well in explaining this rela- 
tionship. However, it seems that some others pulsars possess 
time lag slightly less than the expected value, on average it 
is smaller by an amount of 0.1 in phase. This has probably 
to find its root in the geometry of the polar cap or in retar- 
dation effects not taking into account in this investigation. 

3.6.9 Corrected gamma-ray luminosity 

Finally, knowing the precise geometry from the light-curves 
and time lag, we computed the flux correction factor for each 
of the above mentioned pulsar and summarize our results in 
Tab[T] The "true" efficiency is also given therein. 



4 CONCLUSION 

By combining the striped wind model with a simple prescrip- 
tion for the polar cap geometry and emission, we were able to 
derive the time lag between radio and gamma rays according 
to the high-energy pulse separation, if available, in agree- 
ment with recent Fermi observations from the first gamma- 
ray pulsar catalog. According to our composite model, it 
seems that the observed gamma-ray pulsation is generated 
just outside the light-cylinder and radio pulses come from 
low altitude polar cap locations. 

A further careful inspection of individual pulsar light- 
curves, the correlation between their radio and gamma-ray 
peaks, will permit to severely constrain the geometry, de- 
ducing the obliquity and inclination of line of sight of the 
system. 

The overestimate of time lag by roughly 0.1 in phase for 
many pulsars suggest that another ingredient is still miss- 
ing to fit properly the data. As shown, curved space-time is 
negligible, so we expect magnetic field line bending and/or 
plasma flow within the magnetosphere and other plasma ef- 
fects to give us some clue to this enigma. 
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Figure 8. Gamma-ray (red lines) and radio (blue lines) light-curves of several pulsars fitted with the combined polar cap/striped wind 
model. Intensities are normalized to unity. 



